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1. INTRODUCTION 


À. | BACKGROUND 


As costs for material and energy continue to rise, the 
need for smaller, more efficient heat exchangers continues 
to grow. For those applications involving marine vehicles, 
size and weight limitations dictate the use of small, 
highly-efficient heat exchangers. Reducing the size and 
weight of steam condensers used aboard J. S. Navy ships for 
propulsion and electrical power generation would result in 
lower material costs, and help to aleviate the cramped 
conditions so typical of machinery spaces. Thus, there is a 
strong motivation for continued study in this area. 

The effectiveness of condensers is limited by the water- 
Side, vapor-side, and wall thermal resistances of the 
condenser tubes. Reducing any one of these will contribute 
to increased heat-transfer performance, and smaller physical 
Size of condensers. Methods of enhancing the vapor-side 
coefficient include the use of "roped" tubes, fluted tubes, 
drainage strips attached on the tubes, finned tubes, and 
coatings applied to enhance dropwise condensation. Thi 
thesis concentrates soley on finned tubes. 

Since the late 1940s, externally-finned tubes have been 
used to increase the vapor-side heat-transfer coefficient of 
tubes used in refrigeration systems; but condensers used in 
steam systems, such as Shipboard propulsion plants, continue 
to use smooth tubes. The high surface tension of water, 
which leads to its tendency to flood the area between fins, 
tas resulted ina widely-held belief that such tubes are 
inappropriate for use in steam systems. Recent studies 
[1,2], however, have shown that finned tubes can signifi- 


cantly enhance the heat-transfer rates in such systems. 
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The theoretical treatment of condensation on  finned 
tubes is extremely complex due to the large number of vari- 
ables and physical mechanisms involved. The interaction of 
gravitational and surface-tension forces lead to complex 
three-dimensional flow patterns, which are further dependent 
oğrun Spacing, height, “and thickness. Other variables 
include heat flux, vapor shear, tube diameter, fin Shape and 
fluid properties just to name a few. In view of the above, 
any theoretical models will recuire numerous simplifying 
assumptions, and require complex computer solutions 
involving implicit numerical techniques (3]. To confirm the 
validity of theoretical models, reliable experimental data 
which cover a wide range of relevant parameters must be 
obtained. The availability of such data may then lead to 
the development of fairly simple experimental correlations, 
which could be used in the design stage to predict the heat- 
transfer performance of finned tubes. 

This thesis effort is a continuation of research being 
conducted at the Naval Postgraduate School (NPS) under a 
grant from the National Science Foundation. The basic test 
apparatus used to ccllect experimental data was built by 
Krohn [4]. Graber [5] provided the majority of instrumenta- 
tion, and took preliminary data as the system experienced 
problems with non-condensing gases and partial dropwise 
condensation on copper tubes. Poole [6] made further 
improvements on the apparatus as well as on the instrumenta- 
tion and, most importantly, assured a leak-free apparatus. 
Unfortunately, he did not have sufficient time to produce 
useful data, mostly due to the considerable time spent in 
Systematically locating and fixing leaks, and due to the 
partial dropwise condensation problem that had not been 
solved. Using this system, Georgiadis [7] was finally able 
to obtain complete filmwise condensation on copper tubes. 


He took data on a number o£ finned tubes with fins of 
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rectangular shape, as well as on smooth tubes. The repeat- 
ability of data obtained by Georgiadis demonstrated the 
accuracy of the test apparatus and associated instrumenta- 
tion which was used essentialiy without modification for 
this investigation. 

The overali objectives of the present program at NPS 
includes the testing of: (a) tubes with rectangularly- 
Shaped fins to find the best fin spacing, thickness, and 
height, (b) tubes with various fin shapes to find the best 
geometry that will maximize heat transfer, (c) tubes with 
different  tube-metal thermal conductivity, and (d) the 
effect of vapor shear on finned tubes. Georgiadis [7] 
tested a total of 25 tubes, which were divided into five fin 
spacings (0.5, 1.0, 1.5, 2-0 and ٦۔0‎ 0۱" 9۹69ٔ ٦ 
nesses (0:5; 0.75, 10, sS and 1.5 nm), and two fin 
heights (1.0 and 2.0 mm). At the conclusion of Georgiadis' 
test program, eight new tubes with fin heights of 0.5 and 
1.5 mm, fin spacings of 1.07 1:57 +0 +108 
fin thickness of 1.0 mm, remained to be tested to complete 
the sequence of tubes with rectangularly-shaped fins. The 
testing of these tubes and additional tubes toward objec- 
tives (b), (c), and (d) were the primary goals of this 


thesis effort as listed in the next section. 


B. OBJECTIVES 


The main objectives of this thesis are as follow: 
1. Take data on tubes with rectangularly-skaped fins of 
various fin spacings and fin heights to augment 


previous data [7], 


2. Take data on tubes with fins of different shapes 


(triangular, trapezoidal, parabolic, etc.), 


3. Take data on commercially-available finned tubes, 
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Take data on the “optimum rectangulanly=shaped 


finned tube at different vapor velocities, 


Take data on a stainless steel tube vith "optimun," 


rectangularly-shaped fins, and 


Develop a preliminary correlation based on data for 


tubes with rectangularly-shaped fins. 
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II. THEORETICAL TREATMENT OF CONDENSATION ON HORIZONTAL 


ጨ፦ = Sas s G s Q ee eee a = = — < < re ee ee 


سے سے کے س >= ҹанини АМ.‏ چ — መ-፦=‏ 


A. FILA CONDENSATION 


When filmwise condensation of a vapor takes place on 
smooth tubes, a thin layer of condensate is formed which 
thickens with increasing distance around the perimeter of 
the tube. This condensate layer creates a thermal resis- 
tance, which can limit the heat-transfer performance of the 
tube. This film thickness and its thermal resistance can be 


reduced Гу the use of external, radial fins which, in addi- 


tion to a surface area increase, promote surface-tension 
effects. 
In 1984, Yau et al. [1] neasured the enhancement 


provided by copper finned tubes over smooth tubes for filn- 
wise condensation of stean. Similar experiments by 
Wanniarachchi et al. [2] also in 1984 confirmed that the 
observed enhancements were greater than could be explained 
Ly the increased surface area alone. This additional 
enhancement may be a result of the surface-tension forces 
which act to thin the condensate film. Figure 2.1 schematı- 
cally depicts this phenomenon. 

The effect of surface tension on pressure at the inter- 
face between a liquid and vapor is inversely proportional to 
the radius of curvature of the interface. If a surface has 


two radii or curvature at right angles (i.e., r, and 7د‎ it 


1 
can le shown that: 


1 .2( اما 


Where 


5255 05) ten on cl coldest, dnd 


> 
Ry 
" 


۴۲۰۰۰۰۱۶۲ 50۴6 difference. 


Condensate 









. 


Figure 2.1 Schematic of Condensate Profile at Upper 
Tube Surface. 


| ia 


SC The case of a finned tube, the radius of  curvature 
acound the outer fin perimeter is very iarge compared to the 
radius of curvature around the fin cross-section profile 
TEf yüre 2.1). Referringmback to squation (2.1), it can Be 
seen that the smaller radius of curvature term will domi- 
nate, со surtace-tension effects around the fin perimeter 
may be neglected. 

Because of the convex shape of the condensate film at 
Point A; the pressure within the film at this point is 
greater than the surrounding vapor pressure. Ih a similar 
manner, the pressure at point B is less than the surrounding 


vapor pressure owing to the concave shape of the condensate 
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rın. The relatively flat shape of the condersate film at 
point C leads to a pressure essentially edual to the 


surrounding vapor pressure. These pressures are given by: 


P =p 9 (2220 
ን 
= 6 
ии = (2.3) 
p 
C V (2.4) 
Nene 
P, 7" Surrounding vapor pressure, 
Par Pp Po = liquid pressure at ports AB anc Cee 
Гд, Гә, Гс = radius of curvature of the condensate 


film at ролл Ев рир, али 


At point B, the radius of curvature is small, so the pres- 
sure at point B is less than the pressure at point ር. 
Further, from eguations (oso) an (2.3), the pressure at 
point A aS greater than 1 7٦ In reality, the pres- 
sure gradient within the condensate film varies along the 
height of the fiw due to the continuously varying radi Smo: 
curvature from the fin tip down to the base (3]; however, to 
Simplify the treatment of condensate flow toward the fin 
base, the pressure differences between points A and B and 


points C and B may be simply written as: 


ኪን 
ኪን 


1 ] 
AP = 2 ረ 
ur c E - (25555) 


O 
ከ ከጠ 2486 
ርሾ ፐከ (2.0) 
where 
eevee Ge = prossure difference between points A and B, 


and points C and B. 


De seen, thes pressure differences are positive, 
1 PE in condensate ¿low towardđ»point B. Due to the 
Dome Large mass of Condensate at point 5, gravita- 
tional forces dominate over surface-tensior focces causing a 
O Of condensate around the tube perimeter at the fin 
Lase. In this manner, a condensate run-off channel is 
formed at point B, and the improved drainage of conlensate 
thins the fiim between fins and ou the fin surface. This 
thinning, in turn, reduces the thermal resistance through 
 ን፡ጃ፡ንቢ ን እ እ ሙሮ .፣ ብጋ ። tnus producing an enhancement in adli- 
tion to the gain in surface area. 

It shoull Le stresse1 at this point that the shape of 
the condensate film is highly dependent on the fin geometry. 
Cleamay, then, the fineshape must Le tàken into considera- 
tion in order to maximize the beneficial effects of surface 
tension. 

The above gains may be partially or totally offset, 
however, by tLe tendency of condensate to flood the area 


between fins, especially on the lower part of the tubes. 
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The extent to which the film floods the tube is defined by 
the condensate retention angle (y), as shown in Figure 2.2. 
In the flooded portion of the tube, the relatively thick 
condensate film increases the thermal resistance, thereby 
leading to a degradation in the heat-transfer perfcrmance. 
The condensate retention angle is highly dependent on fin 
spacing, so very closely-spaced fins may lead to completely 
flooded tubes. 

A number of studies have been conducted to investigate 
the performance of finned tubes, and these are presented in 


the next section: 


B. CONDENSATE RETENTION 


In 1946, the first measurements of condensate retention 
were made by Katz et al. [8]. These measurements were made 
under static Conditions Gee. no condensation taking 
place), using water, aniline, acetone, and carbon tetrachlo- 
ride on a number of tubes with different fin densities (276 
to 984 fins/m), and fin heights (122 to 5.-7 mm). It was 
shown that as much as 100% of the tube surface could be 
¿loocec by retained condensate, depending mainly on the 
ratio of surface tension to liguid density and on the fin 
spacing: 

Eight years later, Gregorig [9] recognized that surface- 
tension effects could play an important role on a vertical, 
fluted surface which uses a minimum radius of curvature at 
the flute tip, which gradualiy increases toward the trough. 
As discussed in section A above, this variable radius of 
curvature results in a condensate film pressure that 
decreases toward the trough. In this manner, the condensate 
layer on the convex region of the flute is thinned signifi- 
cantly, thus pushing the condensate into the troughs (Figure 


2.3) where gravitational forces result in enhanced drainage. 
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Figure 2.2 Schematic of Condensate Retention 
on Finned Tubes. 
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Since the Leat-transfer coefficient Gs inversely 7ۃ‎ 
tional to the thickness of the condensate layer, heat 
transfer performance at the trough ~as degraded, = pute ted 
lesser degree than the enhancement experienced by the crest. 
The resulting averaye heat-transfer coefficient, therefore, 


Shows a Significant enhancement. 
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Figure 2.3 Cross Section of the Gregorig Surface 


in 1981, Rudy and Webb [10] made measurements of conden- 
Sate retention angles on finned tubes with three different 
fin densities (748, 1024, 1378 fins/meter) using water, 
R-11, and n-=pentane as Cle ፡፡፡፡፡ሪ ስዕ Later, IN O 
Pudy and Webb [11] developed an analytical model to predict 
the fraction of tube surface that is flooded during conden- 
sation on a horizontal, integral Tin CLs They found that 
the vertical-rise height of condensate on a finned tube vas 
the same as that obtained on a vertically oriented, flat 
finned plate that was obtained by splitting and unrolling an 
identical finned tube. Based on this observation, they used 
capillary equations that predict the liquid rise on a 
vertical U-shaped channel, and assuming negligible vapor 


Shear, they predicted the condensate retention angle to be: 
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(2.7) 


where 
e = fin height, 
er fan thickness, 
5” density of condensate, 
۰و‎ ..eratıon of gravity, 
S - fin spacing, and 
bo = root diameter of tube. 


This result shows that the condensate retention angle 


ቲቲ ቹፒ ሎር ከ with increasing fin density and with 
increasing surface tension-to-density ratio. It should be 
kepta in nuna, however, that the model was based on a 


vertical surface, so it should only be used for angles less 
than about 30 degrees. For angies below this, experimental 
peus و5611‎ the use o£ water, E-11$ R-12, amsmondjas and 
n-pentane were predicted to within 10 percent. It»is month 
mentioning that equation (2.7) kad been derived in 1982, 
تہ‎ ۱11007660 Russian literature, by Rifert | 12). 

Dike sRudy and Webb [10], Owen et al. [13] also recog- 
nized the need to consider condensate retention while 
analyzing finned tubes. In order wo 8618866 ror this; dan 
assumption was made that the condensate retention angle was 
independent of condensation rate, so a static analysis was 


performed. A simple force balance between sucface tension 


2 


and gravitational forces resulted in an >٥ "۹9+ 9 99+٥ 


condensate retention angle as shown below: 


| NE (2.8) 
pe R s De 


where 
0. = overall diameter of fins. 
This equation 15 the sane as eguna tion mp except that 


equation (2.8) is independent of flin thickness (t). 

In 1983, Honda et al. [3] performed a theoretical anal- 
ysis to determine the condensate retention angle. Using an 
iterative numerical scheme, they found the solution to agree 
with equation (2.8) obtained earlier by Owen et al.. Using 
data of their own, Fonda et al. verified a close agreement 
between the predicted and experimental values of the conden- 
sate retention angle. 

Yau et аР. [1] conducted experiments to determine the 
effects of renoving retained condensate from finned tubes by 
instaliing thin metal drainage strips attached edgewise to 
the bottom tube surface. Experiments were conducted where 
the vapor velocity and fin pitch were varied on finned tubes 
with and without drainage strips. They showed a significant 
reduction in the in the condensate retention angle when the 
finned tubes were fitted with drainage strips. Condensate 
retention angles for condensation of steam, ethylene glycol, 
ani E-113 on finned tubes with drainage strips were found to 
fit the empirical relation listed below to within + 10 


percent: 
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V = 251 1 S 1.66 G (2.9) 
рӱ Р, 5 De 


In 1985, Rudy and Webb [14] expanied their 1983 model in 
order to predict the condensate retention angle for hori- 
zontal finned tubes with fins of arbitrary Shape. As 
before, their model vas based on capillary equations that 
predict the amount of liquid rise, and negligible vapor 
shear. The resulting equation for the condensate retention 


angle is given by: 


Eu EC I ہے‎ (2. 10) 
Do Pf R [Cth +s) e - Ap] 
where 
ик ПҸ perimeter of fin cross section, 


fin base thickness, and 


v 
1 


m 
1 


profile area of fin over fin cross section. 


m order to test the validity of equation (2.10), it was 


compared to experimental data for horizontal finned tubes 


with fins of primarily trapezoidal cross section, and fin 
densities ranging from 630 to 1614 fins/m. Steam, R-11, and 
n-pentane were used as the condensing fluids. Equation 


(2.10) was shown to predict the experimental data to within 
a vercent. Note that eguation (2.10) reduces to equation 


(2.8) for fins of rectangular cross section. 


2 


C. THEORETICAL NODES 


Hm 19077 Beatty and Katz (15) performed a number of 
experiments where various refrigerantS were condensed on 
single finned tubes. To predict the film coefficients, they 
Started with the Nusselt equations for condensation on a 
horizontal tube and on a vertical plate, and with the 
respective proportions of horizontal tube area (between 
fins) and vertical fin area, they calculated the total heat 
transfer rate. From this, an average heat-transfer coeffi- 
cient was obtained based on an equivalent tube diameter. 
Their final expression is simply the Nusselt eguation for a 
smooth horizontal tube, but with the tube diameter replaced 
with the equivalent tube diameter. The leading coefficient 


was modified to fit their experimental data as shown below: 


5 3 1/4 1/4 
h = n,sga {kf PF CPF = Py 8 hfe LE (2. 11) 
uf AT DE 
1/4 YA 1/ با‎ 
2o = o 8٦ (2. 12) 
De Aeff | Do مہ۸‎ ин 
Аси Арис лил (2 7٦ 
2 2 
=. (Deccan 
4 De (2. 14) 


where 


surface area of smooth tube, 


d" 
1 


Ext 
[| 


average vapor-side heat-transfer coefficient, 
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ke = thermal conductivity of condensate, 


መ density of vapor, 

heg = specific enthalpy of vaporization, 

He = viscosity of condensate, 

AT = vapor-side temperature drop, 

D. = equivalent tube diameter, 

n = fin efficiency, 

A. = surface area of smooth tube, 

Az = total surface area of finned tube, and 
2. ective area of finned tube. 


Accuracy of this equation was claimed to be better than 11 
percent for a wide variety of nonaqueous fluids. It should 
be noted, however, that this modjel ignores the effect of 
surface-tension forces; SO, ICE IS. valid “only for Love 
surface-tension fluids such as refrigerants. 


Analytical and experimental studies of condensation on 


horizontal tubes with trapezoidally-shaped fins were 
performed by Zozulya, Karkhu, and Borovkov [16,17] some 
years later in the 1970s. Their experiments confirmed the 


need to consider surface-tension forces in addition to 
gravitational forces when developing mathematical models. 
The analytical solutions were based on the following assump- 
tions: 1) the thin condensate film on the fins was treated 
as a laminar boundary layer with a pressure gradient aiong 
the fin profile caused by surface-tension forces; 2) the 
effect of gravitational and inertial forces on the motion of 


the film along the side surfaces of the ‘fins into the 
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Condensate-filled trough was neglected; 3) the motion of 
condensate in the trough area is laminar and produced by 
gravity; and 4) no condensation takes place orn the flooded 
portion of the tube. Differential equations were obtained 
for the height of condensate between fin bases, and tempera- 
ture distribution along the fin height. With the aid ofa 
computer, solutions for these equations were obtained using 


nunerical methods, leading to the expressions below: 


ከ መመር (2. 15) 
Fo AT 
ے7‎ Te 
5 ? cos ? (2. 16) 
2 با‎ - I 4 
Pf ጸ ደ sin T =y 5 
G : 2 + 
12 ug | 70 aera | (2. 19) 
where 
G = condensate flow rate, 
te =E n ۰ ۷۳ 
a - fin semivertex angle, and 
Z = dimensionless depth of condensate 


between fins. 


Although the above eguations appear to be fairly simple, it 


must be kept in mind that rather complex computer solutions 


E 


Seewsecquired to obtain values for 2. Io cheek the validity 
of equation (2. 15), (22,16) m y. experiments were 
performed with four different finned tubes to condense both 
Steam and R-113 vapor. Agreement to within 5 percent was 
reported. 

lun 1977, Nader [18] presented an analytical model to 
predict condensation on a vertical, single fin attached to a 
Horizontal tube. This model differed from those discussed 
above in that the fin temperature was allowed to vary along 
the height of the fin. For the special case where the fin 
temperature is constant, the Nusselt equation was obtained. 
A "condensation efficiency" was then proposed to account for 


a variable fin temperature, and was defined as: 


W- cmm (2. 18) 


3 
pə 7 756 یو‎ 


kr hf AT (2. 19) 
p = t km 
Be (29809) 
where 
C, = constant of proportionality, and 
Kn = thermal conductivity of fins 


With this definition, the rate of condensate formation ard 
the rate of heat transfer (with variable fin temperature) 


were found to be: 


3 


Eh 
3 
ነ. 556 7 E 42.21) 


f he. ] 


1/4 
222 
Q = 1.8856 F F, m AT ( 


where 


rate of condensate formation, and 


ርን 
| 


total heat-transfer rate. 


ርጋ 
II 


Thus, with a table of values for F, and 7 obtained numeri- 
cally, the rate of heat transfer and condensate formation 
could be calculated by a fairly simple expression provided 
the condensation efficiency is greater than 0.8. However, 
for values below 0.8, the accuracy o£ equations (2.21) and 
(2.22) diminishes guite rapidly. 

In 1979, Patankar and Sparrow | 199) trea Geode tiiceca- Ea 
fin attached to a cooled vertical plate or a cylinder as a 
three-dimensional problem. A thin fin was assumed, and 
temperature variations across the thickness were ignored as 
was lengthwise conduction. Additionally, all heat transfer 
was assumed due to condensation only. Compared to the ideal 
fin solution, results obtained by Patankar and Sparrow indi- 
cated that a significant error existed when using the ideal 
fin solution on all but very long fins. Based on their 


results, the following relations were proposed: 
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1/8 1/2 


175: AT 4 kf uf AT p 2 (223) 
06 Һер о ру 了 2 
: ک2‎ | UE 2/22 
2-1 'የቦ Pr e kr e | kr ur AT Et 7- 24) 
- ni hfg g pr” 2 kf e? 
where 
0 2 heat flux, and 


axial coordinate. 


N 
II 


According to Patankar and Sparrow, these equations should be 
valid for most practical applications; however, no experi- 
mental data were available at the time to confirm this. 

In 1980, Rifert [20] analyzed condensation of vapor on 
horizontal finned tubes enhanced by the effect of surface- 
tension forces. In his analysis, he divided the tube into 
flooded and unflooded zones, and solved a two-dimensional 
form of the energy equation for each Zone. The mean heat 
flux was then determined by integrating over each zone and 
the tube length. In cases where condensate is retained in 
more than half of the tube perimeter, Rifert points out that 
a three-dimensional form of the energy equation must be 
used. Solutions to these equations revealed that, in most 
cases, the fin temperature is very nonuniform, so the mean 
integral  heat-transfer coefficient and its correlations 
should not be used. The above-mentioned zone-by-zone anal- 
ysis was, therefore, recommended. 

Based on their study of condensate retention mentioned 


earlier in section B, Rudy and Webb [10] proposed that “the 
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Beatty and Katz [15] model be modified to account for 


condensate retention, as shown below: 





57 57 
S ۔‎ | (2. 25) 


иһеге nos is computed using equation (2.11). This equation 
neglects any heat transfer through the flooded portion of 
the tube. As a result, it was shown to underpredict the 
average heat-transfer coefficient of condensing R-11 by as 
much as 30 percent when a Significant amount of retained 
condensate was present. 

Using the condensate retention angle they had developed 
earlier (see section B), Owen et al. [13] divided a hori- 
zontal finned tube into an upper unflooded section, and a: 
lower flooded section. Unlike Rudy and Yebb [10], however, 
they accounted for heat-transfer through the flooded portion 
of the tube by noting that condensation occurs on both the 
surface of the retained condensate andthe fin tips. An 
effective thermal resistance was obtained over this portion 
of the tube by assuming parallel paths for heat-transfer 
through the fins and retained condensate. The respective 
thermal resistances, which were added in parallel, were then 
added in series with the tube wall thermal resistance. The 


resulting heat-transfer coefficient was then found to be: 


1 0 7 (2. 26) 
一 U m £ 
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where hi is computed using equation (2.11), and 


Dou 
27 NS ۹ 8 ۳۶ 
۱ ገ arr tae a PER INS D پچ مو سی‎ 
i Kaf f 5 Pe “Pe 9 1٣ شب‎ 
k 5. .v€ thermal conductivity 


This equation was found to precict all the available data 
Kalian accuracy of better than 30 percent. 

Rudy and Verb [21] used equation (2.7) a short time 
later while developing an equation for the heat-transfer 
coefficient. The model for this undertaking was based on 


surface-tension-driven radial flow of condensate on the fin 


surface, with a linear pressure gradient along the fin 
height, and gravity-drained flow of condensate in the 
channel between fins. The Nusselt equation for horizontal 


tubes was used for the tube area between fins, while the fin 
surface was treated by an equivalent gravity model developed 
TT TENE et al. (221 and Ruly ( 23) earlrTer. Wsebefore, vapor 
shear and heat transfer through the flooded portion vere 


neglected, yielding a heat-transfer coefficient of: 





3 2 1/4 
m : Pulu M 6” 
Att Do uf AT 
3 1/4 
NE o O SN (n = y) (2. 28) 
Ahr ur €” TA Tn AT T 


Ə 


where 


I = length of tube, 
A, = surface area oi tube زاك‎ TINS ን 


Ac, = fin surface area. 


This expression provided an accuracy of better than 10% for 
condensation of R-11 on short, finely-spaced fins, but the 
accuracy dropped sharply for larger fins spaced further 
apart. This was apparently due to the assumed linear pres- 
sure gradient on the fin surface, and breakdown of the model 
aS gravity forces became dominart. The use of equation 
(2.28) should be restricted to fin densities from 1200 to 
1400 fins/m, ard fin heights of less than 1 mm. 

In 1984, Honda et al. [24] developed a model to predict 
the  heat-transfer coefficient that took surface-tension 
effects and the non-isothermal behavior of fins into 
account. In order to obtain a solution, however, iterative 
finite difference techniques that required a considerable 
amount of computer time were used. Nonetheless, they showed 
that this model was able to predict the average  heat- 
transfer coefficient for most available experimental data 
(which included 11 fluids, and 22 tubes) within 20 percent. 
Predicticns involving the condensation of steam provided the 
largest errors (as much as 407), due in part to the high 
surface tension of water. This reflects the fact that 
surface-tension effects are still not completely accounted 
FOC: Compared to the other models discussed thus far, 
however, the model of Honda et al. shows considerable 
promise in spite of its complexity. 

Th 1965 Adanmek (27) defined a family of condensate 


surface profiles whose curvatures are given by: 


6 


1 _ 6, (E 9 1) ር 
NEL. Y = 
207 | 6i ә 


where 
T rotation angle of normal to condensate film surface, 
S length oí convex surface of fin condensate file, 
S - coordinate along condensate film surface from 


fan tip- 


The parameter E in equation (2.29) characterizes the aspect 
ratio of the fin cross section. A number of E values and 
their corresponding condensate surface profiles are shown in 
Figure 2.3. With the above information, Adamek obtained an 
expression for the heat-transfer coefficient of the fin 


surface as shown below: 


1/4 


R - 2.149 Ff 9 hg, Og 8$ op (E * 1) 2. 30 
Í Sm Uf m A0 (£ + 233 ( ° ) 


These fin profiles were used by Webb et al. -25] and dy 
et al. [26] in 1985 to develop a molel for condensation on 
horizontal finned tubes. This moiel divided the tube into 
flooded, and unflooded regions, and further divided the 
unflooded region into tube (area between fins) and fin 
areas. To predict heat transfer from the unflooded fin 
surfaces, they used the heat-transfer coefficient proposed 
by Adamek. The Adamek fin profile used to approximate their 


trapezoidal fins is shown in Figure 2.4. The Nusselt 
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Figure 2.3 Adamek [27] Condensate Surface Profiles 
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Figure 2-4 Fin Geometry for the Webb et al. Model [22] 
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equation, written in terms of the film Reynolds number, was 
used to predict heat transfer from the unflooded tube 
surface between fins. Writing the Nusselt equation this way 
takes into account the additional condensate thickness that 
results from the drainage of condensate formed on the fin 
surface. For the flooded region, a two-dimensional computer 
code was used to solve for the heat flux into the tube-side 
coolant (452) * which was compared to the limiting case of 
heat flux (qi) assuning zero iin thickness. This ratio 
($ = qp2/31 ) was used along with a linear temperature 
profile across the condensate film to establish a  heat- 
transfer coefficient for this region. The resulting average 


heat-transfer coefficient is given below: 








2 > Ay A£ п – ህ 
ካ = ከ م‎ = |Һ ا‎ + n h t zəbun 5” 
OW 0 | Ne f mər “ua bass (2. 31) 
2 -1/3 
E ህየ 20832 
h. = 1.514 - - Re, ( ) 
5—— 
ə ke 
mE (2. 33) 
Art 
ቤክ ር (Г = П) 
ነ D (2. 34) 


where hẹ is calculated using eguation (2.30), and 


smooth tube area based on the fin diameter, 


= 
Ӱ" 
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Re film Reynolds number, 


4 1 


6 = ratio of tube side 0 7 ا‎ ንክክ 


side heat flux with fins of Zzero thickness: 


A jetailed step-by-step procedure, Cecon oded COC Obtain ma 
a solution of equation (2.31), is provided in Appendix A. 

To allow for heat transfer from the fin tip, and to 
account for decreased heat transfer due to the thick conden- 
Sate film in the trough area, Webb et al. neplaced سک‎ 7۶٤۲ 
S, + t/2- 6,, where 6, is the average condensate film thick- 


m 
ness in the trough area in the unflooded region: 


B = Af )2.35( 


Since the effects due to t/294amd o n 51 0232231634 ٥ 
to cancel each other, Webb et al. reported only a negligible 
change in S,. This model was shown to predict the heat- 
transfer coefficient for R-11 condensing on horizontal 
finned tubes with fin densities of 748, 1024, and 1378 


fins/m within 20 percent. 
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III. DESCRIPTION OF TEST APPARATUS 


A. TEST APPARATUS 


The same test apparatus used by Georgiadis (7] was used 
for this investigation. A schematic of this apparatus is 
502001 In Figure 3.1. Steam was generated in the bciler, 
which consisted of a 304.8 mm (12 in.) Pyrex Jlass section 
containing ten 4000-Watt, 480-Volt Watlow immersion heaters. 
After passing through a 304.8 mm (12 in.) to 152.4 mn 
(6 in.) reducing section, the steam flowed upward through a 
2.48 m (8 ft.) long section of Pyrex glass piping. At this 
point, a 180-degree bend in the piping re-directed the stean 
downward, where after 1.52 m (5 £t.) of Pyrex glass piping, 
it entered the stainless steel test section illustrated by 
751951012 3.2. The test tube was mounted horizontally in the 
test section behind a view port, which allowed visual obser- 
vation of the condensation process. Steam not condensed by 
the test tube continued downward to the auxiliary condenser, 
where it was condensed by two 9.5 mm (3/8 in.) diameter 
water-cooled copper tubes helically coiled to a height of 
ን ከክ (18 in.) The condensate was then returned to the 
boiler by gravity. 

Filtered tap water was used to cool the test tube ona 
once-through basis. This water was first collected ina 
large sump with a capacity of about 0.4 cubic meters (Figure 
3.3), then pumped through a flow meter, and the test tube by 
two centrifugal pumps connected in series. A valve on the 
discharge side of the second pump allowed the velocity of 
water flowing through the test tube to be varied from 0 to 
4.4 m/s (14.4 ft/sec). The auxiliary condenser was cooled 
by a continuous supply of tap water, which was throttled to 


control the internai ressure of the test apparatus. 
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The system used to remove non-condensing gases is shown 
10 ии. The vacuum pump continually drew a sample 
mixture of any possible non-condensing gases and moisture 
from the auxiliary condenser. This mixture was then passed 
through a heat exchanger, where the moisture was condensed 
before it entered the sump, using the same filtered tap 


water described earlier. 


B. INSTRUMENTATION 


The electrical power input to the boiler immersion 
heaters was controlled by a panel-mounted potentiometer. A 
converter with an input voltage of 440 VAC generated a 
Signal which was fed to the data acquisition system to 
calculate the boiler input power. A more detailed descrip- 
tion of the boiler power supply is provided by Poole [6]. 

The internal system pressure was measured manually using 
a U-tube, mercury-in-glass manometer graduated in millime- 
ters. Steam, condensate return, and ambient temperatures 
were measured using calibrated copper-constantan thermocou- 
ples. These thermocouples had an accuracy within + 0.1 K. 
The cooling water temperature rise was measured by two 
Hewlett-Packard (HP) 2804A quartz thermometers, along with a 
10-junction, series-connected copper-constantan thermopile 
as a backup. Throughout this investigation, the quartz 


thermometers and the thermopile agreed to within + 0.03 K. 


C. DATA ACQUISITION SYSTEM 


To monitor the system temperatures and boiler input 
power (using the converter signal), an HP 9826A computer was 
used to control an HP 3497A Data Acguisition System. Raw 
data were processed immediately and, at the same tine, 
stored on diskett for reprocessing at a later time. 
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D. TUBES TESTED 


For this thesis effort, eighteen copper tubes, and a 
stainless steel tube were manufactured. They are presented 
in Table I. The first eight tubes with rectangularly-shaped 
fins were tested in order to complete the sequence of data 
initiated by Georgiadis [7]. This conpleted the systematic 
variation of fin spacing, thickness, and height so that the 
optimum dimensions could be determined. Figure 3.4 shows a 
photograph of the four tubes with fin thickness of 1.0 mn, 
fin height of 0.5 mm, and spacings of 1.0, 1.5, 2.0, and 4.0 
mm (tube numbers 26 to 29). Also shown is the insert used 
to enhance the inside heat-transfer coefficient, which will 
be discussed in more detail in the next chapter. Figure 
3.5(a) shows a cross section of a rectangular fin with a fin 
height of 1.5 mm, a fin spacing of 1.5 mü, and aû TIN ENC 
ness or “nn. In addition to the tubes vith 
rectangularly-shaped fins, six tubes vere tested to study 
the effect of fin shape on the heat-transfer perfcrnance. 
Tubes number 40 and 42 vere made of stainless steel to test 
the effect of fin-metal thermal conductivity on the heat- 


transfer performance. Due to limited machine-shop facili- 
ties available, the "parabolic" fins on two of the tubes 
were not truly parabolic, as shown in Figure 3.5 (b). A 


cross section of a commercially-available finned tube 
produced by the Wolverine Division of United Oil Products is 
also shown jin Figure 3.5(c) for comparison. The smooth 
tubes were tested in order to determine the water-side heat- 
transfer coefficients, and to serve as a comparison for the 
finned tubes. 


E. | VACUUM INTEGRITY 


For any condensation experiment, especially when oper- 


ating under vacuum conditions, a leak-free apparatus is of 
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vital importance. The presence of any leaks would allow the 
admission of non-condensing gases, whose deleterious effects 
would invalidate any data being collected. TO prevent this 
problem, vacuum integrity tests were performed on the test 
apparatus both before and after the data acquisition phase 
of this thesis effort (a time period of about four months). 
The first test revealed a leak rate represented by a pres- 
sure rise of 1.0 mmHg over a period of five days, while the 
pressure rise during the final test was 1.0 mmHg over a 
24-hour period. These tests were performed at an absolute 
pressure of about 85 mmHg. This negligible leak rate, 
together with continuous venting of the test apparatus 
during data runs, eliminated the deleterious effects of 
non-condensing gases on the heat-transfer results reported 


in this thesis. 
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IV. DATA COLLECTION AND REDUCTION 
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A. SYSTEM OPERATION 


In order to ensure filmwise condensation, all tubes were 
treated with a solution containing equal parts of sodium 
hydroxide and ethyl alcohol. The tubes were placed ina 
steam bath and several layers of the solution were applied 
with a tooth brush every 10 minutes. After about one hour, 
the tube surfaces vere completely blackened. The same 
procedure was repeated prior to each data run, but the tube 
was only steamed for ter to fifteen minutes. Following the 
steam bath, the tubes were rinsed with tap water before 
installing them in the test apparatus. Georgiadis [7] 
discusses this procedure in more detail. 

The test apparatus was brought to operating pressure and 
temperature using the procedures provided by Georgiadis | 71. 
Steady-state conditions were assumed to occur once the test 
section steam temperatures, as well as the cooling water 
inlet and outlet temperatures stabilized. Two sets of data 
were then taken for each of the heat fluxes, which were 
determined by the flow rate of cooling water through the 
test tube. Starting at a flow rate of 80 percent (4.44 m/s 
75777 mm O.D. tubes, and 3.92 n/s for 13.7 and 15.80 9m 04D. 
tubes), flow rates were reduced to 70, 60, 45, 35, 26, and 
20 percent, then brought back up to 55. and 80 percent. 
These settings were selected since they provided nearly 
equally-spaced heat flux values. It should be noted that 
continuous adjustments of cooling water flow through the 
auxiliary condenser were required to maintain system pres- 


sure as flow rates through the test tube were changed. 
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Visual observations of condensation taking place on the 
test tube vere made on a regular basis to ensure uniform 
filmwise condensation. A more reliable indication of this, 
however, was the trend of data taken as the cooling water 
Slow rates were increased back up to 80 percent. If the 
final value of the steam-side heat-transfer coefficient at 
80 percent was significantly higher than those at the begin- 
ning of the data run, dropwise condensation was assumed to 
have occurred, and the data were disregarded. All data 
presented in this thesis displayed less than 5% disagreement 
in the steam-side heat-transfer coefficient between initial 


and final data sets. 


B. THE DROPWISE CONDENSATION PROBLEM 


The purpose of this thesis was to take strictly filmwise 
condensation data, using primarily copper tubes. Due to the 
poor wettability of copper with water, especially when even 
minute amounts of contamination are present, there is a 
tendency to condense steam in the partial dropwise mode. 
Dropwise condensation is far more effective than filmwise 
condensation; so, its presence can lead to large errors in 
the data. 

By following procedures set forth by Georgiadis [7], the 


dropwise problem was minimized. During the initial testing 
phase of this thesis, however, a small amount of dropwise 
condensation was encountered. This problem was soon traced 


to the boiler water, which had not been changed for some 
time. By changing the water after every three or four runs, 
the problem was eliminated. 


C. STEAM VELOCITY LIMITATIONS 


A major assumption in the Nusselt theory is that vapor 


velocity does not induce shear forces on the vapor-liquid 
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Miter Lace. Since system pressure in the test apparatus was 
controlled by condensing steam in the auxiliary condenser, a 
reasonable steam velocity was always present. Since stean 
velocity (i.e., steam flow rate) depends on boiler pover 
input, decreasing the power results in lower vapor veloci- 
ties. Under vacuum conditions, (about 85 mmHg) a maxinun 
velocity of about 8 m/s was possible at full boiler power 
(3607 KW). However, at these high velocities, the entire 
apparatus experienced considerable vibration. As discussed 
by Georgiadis [7], a steam velocity of 2 m/s under vacuun 
provided the most stable operating conditions. This corre- 
sponding boiler power (about 9 kW) was low enough to mini- 
mize system vibrations, yet high enough to avoid rapid 
fluctuations in condenser pressure due to the intermittent 
break-up of vapor bubbles in the boiler. To allow a suffi- 
cient amount əf steam to reach the auxiliary condenser when 
operating at atmospheric pressure, the maximum boiler power 
was required. This resulted in a steam velocity of approxi- 


mately 1 m/s. 


D. DATA REDUCTION 


AIT the programs, property®=" functions, "awd caiibpatien 
data used for this thesis effort were essentialy the same as 
those used by Georgiadis [7] and Poole [6]. Because six of 
the tubes (tube numbers 39 thru 44) tested during this 
effort had inside diameters different than any tube tested 
previously, and two (tubes number 40, and 42) were manufac- 
ture fron stainless steel, some minor modifications were 
necessary for the data reduction program (program DRP6) 。 
These modifications included options for different tube 
diameters, tube inserts, thermal conductivity, mixing- 
chamber calibrations, and leading constants for the inside 


heat-transfer coefficients. 
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The mixing-chamber (See Figure 3.2) calibration was 
required to account for the temperature rise that occurred 
across the mixing chamber, due primarily to turbulent 
viscous dissipation. The purpose of the mixing chamber is 
to ensure a uniform temperature at the cooling water 
discharge, so that accurate temperature measurements can be 
made. Since tubes with different inside diameters, 
requiring different inserts, result in different pressure 
drops, calibrations had to be performed for each inside 
diameter. 

A modifiei Wilson Plot program [7] was used to compute 
the leading constants for the Sieder-Tate equation used to 
determine the  water-side heat-transfer coefficient. Data 
taken on smooth tubes were used to determine the leading 
constants, which were then used to reduce data taken on 


similar tubes with identical internal dimensions. 
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V. RESULTS AND DISCUSSION 


A. INTRODUCTION 


During this thesis effort, a number of data runs were 
made using the procedures described in Chapter IV. For both 
vacuum and atmospheric conditions, every tube was tested at 
least three times, on different days, to ensure the data 
were repeatable. Complete filmwise condensation was main- 
tained, and the computed non-condensing gas concentration 
was held below 1.0 percent for all accepted data runs. EG 
was eStimated that an error of + 0.2 K in the steam tempera- 
ture, or an error of + 1.0 mmEg in the system pressure would 
lead to an error of + 1.2 percent in the non-condensing gas 
concentration. This shows that the computed non-condensing 
gas concentration was zero to within the accuracy of temper- 
ature and pressure measurements. As discussed earlier in 
Chapter III, the test apparatus would allow only a negli- 
gible amount of non-condensing gas to be present in the 
apparatus. The non-condensing gas Concentration was 
computed for the system temperature and pressure during 
every data run to ensure no major leaks developed. A 
summary of tubes tested between Georgiadis 77) andü this 
thesis effort and the resulting enhancements are provided in 
Table II. 


B.  WATER-SIDE HEAT-TRANSFER COEFFICIENTS 


The water-side heat-transfer coefficients, enhanced by 


the use of an insert, were obtained by Georgiadis [7] for 
the tubes with an inside diameter of 12.7 mn. Both a 
"direct method," and a "modified Wilson method" were used, 


and are discussed beiow. 
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1. The Direct Method 


For the "direct method," a smooth tube was manufac- 
tured with six wall thermocouples inserted around its perim- 
eter. During data runs under vacuum and atmospheric 
conditions, readings from these six thermocouples were 
obtained and averaged, in order to determine an average tube 
wall temperature for each heat flux. It should be noted 
that the wall temperature was guite non-uniform around the 
tube perimeter. Georgiadis [7] reported up to an 18-K drop 
in the wall temperature from the top to the bottom of the 
tube. Evaluation f the uncertainty resulting from this 
temperature variation was considered beyond the scope of the 
investigation. The average wall temperature was then used 
along with the cooling water inlet and outlet temperatures 
to determine the log-mean-temperature difference  (LMTD). 
With values for heat flux and LMTD known, the value of the 
water-side heat-transfer coefficient was obtained. This, in 
turn, was used to determine the leading coefficient for the 


Sieder-Tate-type equation shown below: 


Orr 
al 
69. 5. n, 1/3 Ë 十 了 (5. 1) 


Ну 


A more detailed step-wise solution procedure is given by 
Georgiadis [7]. For th2 previously-mentioned tubes of 
15.8 mm and 19.05 mm O.D., Georgiadis determined the coeffi- 
cient C TG Le 0 06 ) ٣ The constant B = 26.6 15 an addi- 
tional parameter found to improve the fit of equation (5.1). 
Note that C = 0.0635 is about 2.5 times greater than the 
well-known Sieder-Tate coefficient of 0.027 [28] for long 
tubes with smooth inner walls. This is due primarily to the 
insert used to enhance the water-side heat-transfer 


coert1cienes 
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D The Modified Wilson Method 
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À complete description of this method is provided in 
Appendix B. Briefly, this method uses a Sieder-Tate-type 
equation for the water-side heat-transfer coefficient, and 
an equation (discussed below) for the steam-side  heat- 
transfer coefficient. Both equations contain leading coef- 
ficients that must be determined iteratively while being 
used in a linear equation to generate the Wilson plot. The 
iterations were performed on the HP 9826A computer mentioned 
previously in Chapter III, using program "WILSON6." 

The "modified Wilson method" used during this inves- 
tigation consisted of one minor modification from that of 
Georgiadis. For this thesis, a Fujii-type [29] equation was 
used instead əf the Nusselt-type equation for the steam-side 
coefficient. Ihis modification was made to account for the 
small steam velocity (about 1 m/s) that was present during 
the runs. For this reason, use of the Fujii-type equation 


appears to be more accurate than the use of the Nusselt-type 


equation. This modification resulted in a slightly higher 
value (up to 3 percent) for the leading coefficient C. The 
program "WILSON6" described earlier, allows an option for 


selecting eithér the Fujii-type or Nusselt-type equation for 
the steam-side coefficient. 

Using a Nusselt-type equation for the steam-side 
coefficient, Georgiadis (7) found the leading coefficient 
(some eqwation (2. 11) C to be 0.071, with the B value set 
ecual to zero. This C value is about 10 percent higher than 
the value (C = 0.0635) obtained by the "direct method," but 
Since the "direct method" is generally felt to be more reli- 
able, the values of 0.0635 and 26.4 were used for the 
constants C and B respectively for this thesis effort with 


tubes of 12.7 mm inside diameter (see Table I). 
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Due to the thin (i.e., 0.5 mm thickness) tube walls 
for tubes 39 thru 42, it was not possible to manufacture an 
instrumented tube. For this reason, the "modified Wilson 
method," together with the Fujii-type equation for the 
steam-side coefficient, was used to determine the leading 
coefficient for equation (5.1) for both the copper tube 
(tube number 41), and the stainless steel tube (tube number 
42). For the copper tube, a value of C = 0.0756 was 
obtained, while for the stainless steel tube, a value of 
C = 0.0688 was obtained. Although the inner tube diameters 
for both tubes are the same (i.e., 13.5 mm), the much higher 
thermal conductivity of copper, which results in circunfer- 
ential heat conduction through the tube wall, leads to its 
higher value for C. A summary of the leading coefficients 
used in equation (5.1) for all the tubes is presented in 
Table III. 


C. REPEATABILITY OF DATA 


Since the reliability of data taken was of vital impor- 
tance, all data runs were repeated, as Georgiadis [7] did, 
at least three times on different days. The computed stean- 
sile coefficients for similar conditions (i.e., same tube 
and the same operating conditions) from different days 
agreed to within + 5 percent. Additionally, data runs were 
performed on two finned tubes (tubes number 6 and 17), under 
Similar conditions to verify the repeatability with data 
taken by Georgiadis. Figure 5.1 shows that the experimental 
steam-side heat-transfer coefficients of Georgiadis and 
those obtained during this investigation agree to within + 5 
percent. The curves shown in this figure (and all other 
figures of a similar nature that follow) are the least- 


squares-fit curves according to the following equation: 
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5 ATP (5.2) 


where a and b are experimentally determined constants. 
These curves were selected instead of those based on the 
steam-side coefficient versus heat flux since the uncer- 
tainty band on AT is more uniform than the uncertainty band 
based on the steam-side coefficient. À brief discussion of 
the uncertainty bands is presented in Appendix D, along with 
a listing of the uncertainty analysis program  "UNA6" and a 


few sample runs. 


D EFFECTS OF FIN SPACING AND FIN HEIGHT ON PERFORMANCE 


Data were taken on eight tubes (see Table II) with fins 
of rectangular shape (tubes number 26 thru 33) under vacuum 
(approximately 85 mmHg) and atmospheric conditions. Four 
tubes had a fin height of nn ands fin Spacings. of 
ii 1.5 2.0, and 4.0 mn, while the remaining four tubes 
had a fin height of 1.5 mm, and the same sequence of fin 
spacings. Georgiadis [7] reportei data on rectangularly- 
shaped fins with the same fin spacings, but with fin heights 
of 1.0 mm, and 2.0 mm. A complete listing of tubes tested 
by Georgiadis is also presented in Table II (tubes 171 thru 
29) 

Data reductions were performed on the EP 9826A computer, 
using a program named "DRP6." This program was essentially 
the same one used by Georgiadis ("DRP4Y"), with modifications 
to allow for variations in tube diameters, and tube 
material. Data obtained for all tubes tested during this 
thesis effort (tubes number 26 thru 44) are presented in 


Appendix C. 
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This section presents results showing the variations 
of the steam-side heat-transfer coefficient with heat flux 
having fin spacing aS a parameter. Figures 5.2 Vand 7 
present data for tubes with a fin height of 0.5 mm under 
vacuum and at atmospheric pressure respectively. The smooth 
tube data and the theoretical prediction of Nusselt are are 
included for comparison. In a Similar way, Figures 5.4 and 
5.5 present data fox tubes with a fin height of 1.5 mnm. 
These figures show that the best heat-transfer performance 
occurs for the tubes with a fin spacing of 2.0 mn, while the 
worst performance occurs for the tubes with a fin spacing of 
4.0 mm. Note that these trends are the Sane Езел ии 
vacuum and at atmospheric conditions. 

The comparison of finned tubes is made through the 
enhancement ratio. This ratio is defined as the steam-side 
heat-transfer coefficient of a finned tube to that of the 
smooth tube (same diameter as the finned tube root diameter) 
at the same heat flux. As was Shown by Georgiadis [7], the 
uncertainty in calculating the steam-side heat-transfer 
coefficient increases with decreasing heat flux; so, heat 
flux values of 0.35 and 1.0 MW/m? were chosen for vacuum and 
atmospheric conditions respectively. Cross plots of the 
enhancement ratio versus fin spacing are shown in Figures 
5.6 and 5.7 Tor fin heights of 075, 0 777771 MU 
(curves representing e - 1.0 and 2.0 mm are from Georgiadis 
[TA Enhancements as high as 4.8 and 6.4 can be seen for 
vacuum and atmospheric conditions respectively. Once again, 
fin spacings of approximately 2.0 mm for the tubes with “in 
heights of 0.5 and 1.5 mm are shown to provide the optimum 
performance. Data provided by Georgiadis 7] for fin 
heights of 1.0 and 2.0 mm, however, demonstrate an optimum 


heat-transfer performance for fin spacings of 1.5 mn. As a 
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= (111, there does not appear to be a clean oend forthe 
optimum fin spacing as a function of fin height. 

On one hand, the difference in performance between 
Bül — 1.5 and 2.0 mm finned tubes, as well as the s = 0.5 
and 1.0 mm finned tubes is less than the experimental uncer- 
tainty (+4 15 percent), so that the optimum spacing may be 
somewhere between 1.5 and 2.0 mm. On the other hand, each 
data point in these figures has been repeated within 5 
percent (1.e., computed steam-side heat-transfer coeffi- 
cient) on different days. Also, as discussed earlier, the 
data taken on tubes 6 and 17 by Georgiadis were repeated 
within 3 percent during this investigation. Therefore, the 
slightly disagreeing trends shown by Georgiadis and by this 
investigator (see Figures 5.6 and 5.7) may be justified by 
the very complex nature of the condensation phenomenon on 
finned tubes. 

The presence of an optimum fin spacing, however, is 
easy to understand. The 1.0-mm fin spacing, which provides 
the largest increase in surface area, performed worse than 
the 2.0-mm fin spacing due to the relatively large conden- 
sate retention angle with the 1.0 mm fin spacing. Visual 
observations showed that the condensate retention angle was 
around 110 degrees for the tube with s = 1.0 mm, while it 
was about 40 degrees for the tube with s = 2.0 mn. The 
additional thermal resistance induced by this layer of 
condensate overpowers the benefit gained from increased 
surface area, so the heat-transfer performance is reduced. 
As fin spacing increases, the condensate retention angle 
decreases, as does the thickness of the condensate film 
between fins on the upper portion of the tube. Beyond fin 
Spacings of 2.0 mm, the decrease in surface area tends to 
overshadow the other effects, resulting in decreased heat- 
transfer performance compared to the tube with the "optimum" 


spacing. 


Ta 


As mentioned earlier, fin heights of 0.5 and 1.5 mn 
were examined during this thesis effort, and data for the 
tubes with fin heights of 1.0 and 2.0 mm were available from 
Georgiadis (7). Table II shows that the general trend is an 
increase of enhancement ratio as fin height increases, which 
results from the increase of surface area provided by the 
higher fins. 

To study the enhancement obtained beyond the 
increase in surface area due to finning, two additional 
columns are provided in Table Il. These columns represent 
enhancement ratio/area ratio (Eo/Ar) at both vacuum and 
atmospheric conditions. As can be seen, Eo/Ar values range 
from 0.86 to 2.23 under vacuum coniitions, and from 1.43 to 
3.20 at atmospheric pressure. Figures 5.8 and 5.9 present 


plots of of Eo/Ar as a function of fin spacing under vacuum 


conditions, and at atmospheric pressure respectively. The 
relatively small values of Eo/Ar - 0.86 (vacuum) and 
Eo/Ar = 1.43 (atmospheric) for tube number 15 are due to tne 


small fin spacing which results in a completely flooded 
tube. On the other hand, the values of Eo/Ar = 2.23 
(vacuum) for tube number 14, and Eo/Ar = 3.20 (atmospheric) 
for tube number 28, Show that significant enhancements 
beyond the area ratio are obtainable when the amount of 
condensate flooding is kept small. It appears that surface- 
tension forces are mostly responsible for this enhancement 
beyond the area ratio. 

As discussed in Chapter II, the condensate thinning 
provided by surface-tension effects is related to the pres- 
Sure gradient from the fin tip to the fin root. For short 
fins, this pressure gradient is relatively larger than for 
high fins, resulting in improved thinning of the condensate 


mum. On the other hard, as fin height increases, surface 
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area increases and the condensate retention angle decreases. 
Both of these effects lead to increased heat-transfer 
performance, but when combined with the poorer thinning 
effect, the Eo/Ar ratio decreases somewhat for the larger 


fin heights. 


E. EFFECT OF FIN GECMETRY ON PERFORMANCE 


In order to study the effect of fin geometry, data were 
taken on tubes with fins of triangular, trapezoidal, and 
parabolic shapes (tubes number 34, 35, 37, 38), as well as 
on a spirally-finped tube with threads of triangular shape 
(tube number 36). Fin dimensions (i.€e., fin spacing, 
height, and thickness) for these tubes were chosen so that 
the resulting total surface areas would be approximately 
equal to that of the "optimum" tube with rectangularly- 
shaped fins and a fin height of 1.0 mm (tube number 6). Bn 
this manner, the effects of increased surface area were 
eliminated, allowing a comparison between tubes where fin 
geometry was the only variable. A commercially-available 
"Wolverine" tube (tube number 43), whose fin geometry is 
shown in Figure 3.5(c), was also tested. Dimensions for the 


above fins are given in Table II. 
1. Effect of Fin Shape 


This section examines the effect that fin shape has 
on the heat-transfer performance of finned tubes. The 
performance of tubes number 34, 35, 36, 37, and 43 under 
vacuum conditions is shown in Figure 5.10, whiie Figure 5.11 
Gepicts their performance at atmospheric pressure. The 
dashed curve in these figures represents a least-squares fit 
for data taken on the "optimum" tube with rectangularly- 
Shaped fins and a fin height of 1.0 mm (tube number 6). In 


both cases, the "Wolverine" tube (tube number 43) gave the 
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worst performance. This was due primarily to the small fin 
Spacing (0.7 mm) that resulted in a completely flooded tube. 
Based on the results of Georgiadis [7], the small fin thick- 
ness (0.3 mm) of the Wolverine tube may also be responsible 
for its poor performance. Georgiadis showed that for a fin 
spacing of 1.0 mm the tubes with a fin thickness of 0.5 mm 
(tube number 16) performed 16 percent and 6 percent poorer, 
under vacuum and at atmospheric pressure respectively, than 
the tube with a fin thickness of 1.0 mm (tube number 5). 
Based on this observation, the small fin thickness is also 
responsible for the poorer performance of the "Wolverine" 
tube in comparison, for example, to the tube with fins of 
parabolic shape (tube number 37), which was also completely 
flooded by condensate. Therefore, the poor performance of 
this tube is probably due more to increased condensate 
Flooding and smaller fin thickness than to fin shape. On 
the other hand, the best performance was observed for the 
spirally-finned tube (tube number 36) under vacuum condi- 
tions, andthe tube with parabolic fins, where fin Spacing 
was 0.7 mm (tube number 37), at atmospheric pressure. The 
fact that tube number 37 performed well at atmospheric pres- 
sure, but poorer than the spirally-finned tube under vacuum 
conditions was somewhat unexpected. It must be kept in mind 
that these data runs were repeated at least three times on 
different days, thus confirming this trend. One possible 
explanation centers on variable fluid properties, as 
discussed below. 

Under vacuum conditions, the condensate retention 
angle is significantly smaller for tube number 36 (about 30 
degrees) than for tube number 37 (completely flooded), and 
the saturation temperature is about 50 K lower. Therefore, 
the relatively high heat-transfer coefficient of the 
Spirally-finned tube is due mainly to its small condensate 
retention angle. When the condensing temperature is raised, 
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both surface tension and viscosity decrease. For example, 
the surface tension decreases by about 10 percent, and the 
viscosity decreases by about 45 percent when the temperature 
is raised from 40 °C to 90 °C (assumed approximate film 
temperatures). Decreased viscosity should result in higher 
heat-transfer coefficients at atmospheric pressure than 
under vacuum Conditions for a given tube. Further, it 7 
be possible that the effect of viscosity is dependent on the 
fin shape. Unfortunately, no equations exist to support 
this hypothesis. The effect of surface tension, on the 
other hand, is fairly well understood; the condensate reten- 
tion angle decreases with decreasing surface  tensicn. 
Therefore, the reduction in retention angle (from vacuum to 
atmospheric pressure) for the tube with "parabolic" fins is 
greater than the reduction for the spirally-finned tube, 
Since the latter was already small under vacuum conditions. 
Thus, the reduced retention angle on the tube with "para- 
bolic" fins provides an increased heat-transfer performance, 
Surpassing the performance of the spirally-finned tube. 

The performance of the remaining two tubes (tubes 
number 34 and 35) is the same under vacuum conditions, but 
at atmospheric pressure, the tube with triangularly-shaped 
fins (tube number 34) outperformed the tube with trapezoidai 
fins (tube number 35) by about 10 percent. An argument 


Similar to the one above may explain this behavior. 
2. Effects on Enhancement 


Referring to Table II” it can be seen that all of 
the above tubes provided significant heat-transfer enhance- 
ments over smooth-tube values. Note that the enhancement 
obtained under vacuum conditions for tube number 36 was the 
same as obtained for the "optimum," rectangularly-shaped 
finned tube with the same fin height (tube number 6). Under 


atmospheric conditions, however, the tube with fins of 
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parabolic shape (tube number 37) performed significantly 
better than any of the tubes with rectangular fins of the 
Same height (1.0 mm). 

Like the tubes with rectangular fins, most of these 
tubes displayed enhancements above the area ratio, as shown 
by the two Eo/Ar columns in Table II. These ratios range 
from 0.94 to 2.95 under vacuum conditions, and from 1.52 to 
4.13 at atmospheric pressure. Like tube number 15, the 
small values of 0.94 and 1.52 for the "Wolverine" tube 
resulted from the large amount of retained condensate, and 
by the very thin fins (about 0.3 mm) as discussed earlier. 
As was the case for rectangular fins, the highest Eo/Ar of 
2-95 and 4.13 were shown for the tube (tube number 36) with 
a relatively small area ratio of 1.37, indicating once again 


the importance of surface-tension effects. 


3. The Performance of "Parabolic" Fins 


In order to study the parabolic fin more closely, a 
tube with parabolic fins was manufactured (tube number 38) 
with the same fin height, spacing, and thickness as tube 
number 17 (tested by Georgiadis [7] and this investigator). 
The performance of these two tubes 1S compared in Figure 
5.12 under vacuum conditions, and in Figure 5. 13 at atmos- 
pheric pressure. Clearly, the parabolic fins provide a 
greater enhancement than the rectangular fins. The reason 
for this, as mentioned earlier, is probably the continuous 
change of radius of curvature (increasing from the trough 
area to fin tip) for the parabolic’ shape. The resulting 


pressure gradient along the fin height is more favorable 


(than on rectangular fins), which, in turn, leads to 
improved drainage ‘from the fin surface. This effect was 
seen most clearly on tube number 37. Another observation 


(that contributes to improved heat transfer) made on the 
tube with "parabolic" fins is that, although  conpietely 
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flooded, the fin tips (about half the fin height) were 
visible along the majority of the fin circumference (except 
for about the bottom 30 degrees). 

AS discussed in section ES dS tube number 37 
performed poorer than the spirally finned tube (tube number 
36) under vacuum conditions. The reasons for this are not 
fully understood; so, further tests of these fin geometries 


(1-e., spiral and parabolic) should be performed. 


me HETFECİ OF FIN-METAL THERMAL CONDUCTIVITY OW STEAM-SIDE 
COEFFICIENT 


To investigate the effect of fin-metal thermal conduc- 
tivity on the heat-transfer performance, a stainless steel 
tube (tube number 400) was manufactured with the sane fin 
dimensions as the "optimum" copper tube (tube number 6) with 
rectangular fins, and a fin height of 1.0 mm. Due to the 
low thermal conductivity of stainless steel, it was neces- 
sary to reduce the tube wall thickness to 0.5 mm to prevent 
an excessive wall thermal resistance. A copper tube (tube 
number 39) with the same dimensions as the stainless steel 
tube was manufactured for comparison. Also, two smooth 
tubes having an outside diameter equal to the root diameter 
of the finned tubes (tubes number 41 and 42) were manuíac- 
tured; the first tube was made of copper, and the other was 
made of stainless steel. The results for data runs taken 
under vacuum conditions, and at atmospheric pressure are 
Shown by Figures 5.14 and 5.15 respectively. 

As expected, the copper finned tube provided a signifi- 
cant enhancement, while the stainless steel finned tube 
actually reduced the heat-transfer performance slightly in 
comparison to the smooth tube (tube number 42). The thermal 
conductivity of stainless steel is much lower than that of 


copper  (i.e., 15 W/m.K compared to 385 W/m.K), so the 
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thermal resistance of the stainless steel fins is signifi- 
cantiy higher, resulting ina poorer fin efficiency. This 
is combined with the deleterious effects of retained conden- 
Sate (introduced by the use of fins) so that the above 
factors dominate over the beneficial thinning of the conden- 
sate film between fins. The net effect, then, isa slight 
degradation of the heat-transfer performance compared toa 


smooth tube. 


G. EFFECT OF STEAM VELOCITY 


As discussed in Chapter II, most models of film ccnden- 
Sation on finned tubes have neglected the effect of vapor 
shear. For this reason, a decision was made to perform runs 
for the "optimum" tube (tube number 6) with rectangular 
fins, and fin height of 1.0 mm, at nominal steam velocities 
of 2 4, 06, amd 09/75: Steam velocity through the test 
section was controlled by the boiler input power, and since 
a Steam velocity of 1 m/s was the maximum obtainable at 
atnospheric pressure, the data runs were performed under 
vacuum conditions only. The results of data taken for the 
above-mentioned steam velocities are shown in Figure 5. 16. 

The trend shown in Figure 5.16 indicates that increased 
steam velocity (i.e., increased vapor shear) leads to 
enhancements of up to 10 percent as steam velocity increases 
from 2 to 8 R/S This is to be expected since the shear 
stress at the liquid-vapor interface tends to pull the 
condensate toward the bottom of the tube, thereby thinning 
the condensate film at the top part of the tube. 

The above-mentioned data were compared to data taken by 
Yau et al. [1] for steam velocities of 0.5, 0.7, and 
1.1 m/s at atmospheric pressure. They reported enhancements 
of up to 40 percent for the steam velocity of 1.1 u/s 
compared to the 0.5 m/s velocity. The rather large 
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difference in enhancements reported by Yau et al. compared 
to those of this investigation was not expected. It should 
be kept in mind, however, that the data of Yau et al. were 
collected at atmospheric pressure (compared to vacuum condi- 
tions for this investigation), where the increased satura- 
tion pressure leads to a lower condensate viscosity (about 
45 percent as discussed in section E.1). This, in turm, 
allows condensate to flow more freely so that the effect of 
shear stresses acting on the liquid-vapor interface would be 
more pronounced. Another factor to be considered is the 
slightly larger fin spacing used by Yau et al. (2.0 mm 
compared to 1.5 mm for this investigation). If the fin 
spacing 1s very small, the vapor shear may not be "felt" on 
the sides of tha fins or on the fin bases. This is due to 
the fact that the steam flows around the tube, thus leaving 
nearly stagnant steam between fins. As the fin spacing 
increases though, the fin sides and bases will be more 
"accessible" to steam, therefore, giving rise to the benefi- 
cial effects of vapor shear. Thus, the larger enhancements 
of Yau et al. [1] compared to this investigation may be 
explained by the two phenomenon just described. 

An attempt was made to correlate the data using a 


rujii-tvpe [29] equation given by: 


Nu = c Fd 


Ferp j (5.9) 


where 


E ۳ 


. (5.4) 


90 


Be z Pf Vs Do 
tp с (512) 


F = g Do علا‎ һер 


Vee kp AT (5.6) 
7 stean velocity 
h g - Steam-side heat-transfer coefficient based 


on actual surface area (including fins) 


Retp - two-phase Reynolds number 


C, 4 = constants to ke determined experimentally 


E constants c and d used in the Fujii equation [29], which 
was developed to account for vapor shear on smooth tubes, 
were 0.96 and 0.2 respectively. Values of c= 1.302 and 


d = 0.236, however, were found to more closely fit the data 


of this investigation after a least-Squares E1 የ. 
perforned. A comparison of the experimental data to this 
cugve is shown in Figure 5.17. Although a fairly close 


agreement has been obtained, the orientation of data clus- 
ters does not match the slope of the least-squares line 
generated by globally fitting ali the data. It is possible 
that the errors in the experimental constants used in the 
Sieder-Tate-type equation are partly responsible for this 
questionable trend.  Georgiadis [7] showed that using larger 
values (from the "modified Wilson method") of the 
Sieder-Tate coefficient (C) tends to rotate the data clus- 
ters to more closely match the experimental correlation. 
Further, a Fujii-type correlation may be inappropriate for 


finned tubes due mainly to condensate retention that occurs. 
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Note that the slope of the line predicted by equation 
1565) zor the experimental data does not match the sloçe o£ 
the Fujii equation mentioned earlier, but instead, nearly 


matches the slope of the Nusselt equation for smooth tubes 


given by: 

3 2 1/4 

h = 0.728 | 5f PF & hfg 

hf Do AT (5779 

which can be rewritten as: 
1/4 
775257 i | 
Re. 7 (5.8) 


etp 


In this manner, the effect of vapor shear, at least on the 
tube tested (tube number 6), is very small compared to 
smooth tubes [7]. In fact, using the experimentally-found 
d-value in equation (5.3), the steam-side coefficient and 


steam velocity are related as follows: 


0.028 
)5.9( ا و 


Once again, this observation shows that the effect of steam 
velocity on the steam-side heat-transfer coefficient is very 
snall. 

A vertical shift of the least-squares-fit compared to 
the Nusselt theory results from the insreased c-value (1. 302 
6ت‎ 0 0 307600-06 0.728) of equüuatron (5:8y- Tn fact, the ratio of 
the experimental c-value to the coefficient in the Nusselt 
equation (i.e.,  1.302/0.728 - 1.788) represents the heat- 


transfer enhancement in excess of the area increase due to 
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finning. Therefore, the c-value depends fin geometry, and 
surface tension of the condensate etc.. 

In view of the above discussion, further attempts to 
molify the Nusselt equation to include vapor shear, surface- 
tension effects, and fin geometry etc. may leai to an accu- 
rate correlation TOL the steam-side heat-transfer 


coefficient for finned tubes. 


H. EFFECT OF INTERNAL AND EXTERNAL ENHANCEMENTS ON THE 
OVERALL COEFFICIENT 


Up to this point, all discussions have delt with 
enhancements to the steam-side heat-transfer coefficiert. 
Even with large steam-side enhancements, however, a large 
thermal resistance through the tube wall, or on the water 
side may dominate, and completely negate the steam-side 
enhancement. For this reason, a decision was made to test 
the "optimum" tube with rectangular fins (fin height of 
1.0 nm - tube number 6), and the corresponding smooth tube 
(tube number 1), with and without the insert. The results 
of these data runs, taken under vacuum conditions, are shown 
in Figure 5.18. At a water velocity of about 3 m/s, use of 
the insert in the smooth tube resulted in a 40 percent 
increase of the overall coefficient, as did the use of 
external fins without the insert. By combining the use of 
fins and the insert, an enhancement of nearly 155 percent 
was attained. It must be noted, however, that the insert 
used during these experiments is not suited for use in 
condensers due to the very high pressure drop it creates. 
Instead, Webb et al. [30] showed that the use of internal 
flutes, fins, or circumferential ribs may be considered to 
increase the water-side coefficient by a factor of up to 2.0 
(compared to 2.5 obtained by the insert used during this 


investigation). As a result, it appears that the use of 
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external fins  togather with internal enhancements, could 
increase the overall heat-transfer coefficient by approxi- 


mately 100 percent. 


I. DEVELOPMENT OF AN EXPERIMENTAL CORRELATION 


Using the 1985 model of Webb et al. [25] (discussed in 
Chapter II) as a starting point, attempts were made to 
correlate the experimental data obtained for tubes with 
rectangularly-shaped fins. Time limitations prevented an 
in-depth effort, but a fev simple modifications were made. 
Before discussing these modifications, the general scheme of 
the Webb et al. model will be discussed. A more detailed 
step-by-step discussion is provided in Appendix A. 

First, a value for the parameter 5 must be determined 
for the fin geometry in question. With this value Known, 
the heat-transfer coefficient for the unflooded fin surfaces 
can be obtained using equation (2.30). Values for the heat- 
transfer coefficients of the unflooded tube area between 
fins, and for the flooded region are determined next using 
equations (2.32) and (2.33) respectively. The above- 
mentioned heat-transfer coefficients are then substituted 
into equation (2.31) to obtain the average heat-transfer 
coefficient for the finned tube in question. 

A careful review of the Webb et al. model revealed two 
deficiencies worth noting. First, in their two-dimensional 
analysis of the flooded tube region, the heat-transfer coef- 
ficient for the flooded root area between fins was applied 
to the fin tips. Figure 2.4 shows that a constant,  rela- 
tively thick condensate film is assumed which, when applied 
to the fin tips, results in an underprediction of the heat 
transfer. Use of the Nusselt eguation for smooth tubes 
would be more appropriate for fin tips in this region. 


Additionally, the effects of convection in this region are 
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ignored; however, to include these effects would make a 
solution extremely complex. 

The second deficiency, probably valid only for 
rectangularly-shaped fins, involves the assumed condensate 
film shape between fins on the unflooded tube surface. 
Figure 2.4 shows a depression in the film next to the fia 
base, with an increasing condensate thickness away from it. 
Condensate film shapes observed during this investigation 
continued to decrease away from the fin base as shown by 
Figure 2.1, probably due to the thinning effect of surface- 
tension forces. Based on the assumption of a relatively 
thick condensate film, predictions using the Webb et al. 
model will underpredict heat transfer from this part of the 
tube. 

To use the Webb et al. model for rectangular fins, a 
value of -0.9 was chosen for zg. A value of -1.0, which 
corresponds to a straight line  (i.e., straight side of the 
rectangular fin) would seem more appropriate; however, this 
value would lead to zero heat transfer from the fin sides 
using equation (A.1). Webb et al. used values around -0.7 
to -0.8 for their slightly trapezoidal fins; so the value of 
-0.9 was assumed to be reasonable. Using this value along 
with a » value of 5.0 (used by Webb et al.) resulted in 
significant underprediction of the experimental data 
obtained by this investigator and by Georgiadis [7], as 
Shown in Figure 5.19. 

In an attempt to improve the predictions, values for the 
ф ratios were calculated using the paralled heat flow paths 
proposed by Owen et al. [13] discussed earlier in 
Chapter II. These ç values are dependent on fin spacing, 
and should decrease with increasing fin spacing. This trend 
was observed, and for the set of tubes with both fir height 
05 Os 2.0, and 


4.0 mm = tubes number 4 to 38), Tesulted in $ values that 


and fin thickness of 1.0 nm (s 
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banged from 1.7 (s = 0.5 mm) to 1.2 (s = 4.0 ma). Using the 
Vebb et al. model once again, but with a constant average 


value of 1.5, gave the results shown in Figure 5.20 for the 


above-mentioned tubes (tubes number 4 to 8). This Figure 
also shows a plot of the Webb et al. model using the vari- 
able (ieee. / 061.2) whachenade, almost»anosdaf ference. 


Once again, the experimental data, also shown in Figure 5.20 
were significantly underpredicted. 

In order to correct for the assumed thick  condensate 
layer between fins, the work of Fujii et al. (31) vas used. 
Their paper delt with smooth tube enhancements from coiling 
a wire around the tube circumference, and provided an equa- 
tion to account for the thinning effect between wire wraps 


as shown below: 


1/4 
Nue s S & (1 + A) (5- 10) 
Мис (s t t) 3 
where 
| ... 4 0 Da (5.11) 
— 
Pf م‎ 5 Pg 
7” 
25 -2 
r_ = 0.03 n 2 
1 = i Y | i 
Pre diameter (taken as the fin height 


for use in the Webb et al. model) 


Adding this correction factor to the heat-transfer coeffi- 
cient for the unflooded root area between fins made a 


significant improvement for the group of tubes) mentioned 
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earlier (tubes number 4 to 8) as shown by Figure 5-20. 
Nonetheless, experimental data were still underpredicted. 

In summary, it should be kept in mind that this model 
was not developed for rectangular fins. The most signifi- 
cant discrepancy, however, appears to be on the prediction 
of heat transfer through the flooded area. This model 
predicts only minimal heat transfer through the floodec 
portion of the tube, whereas Wanniarachchi et ai. (27 
reported enhancements at least equal to the area ratio for 
conpletely flooded tubes. Additionally, there may be other 
phenomena which have not been taken into account (such as 
intense convection in the condensate) which, if included, 
may inprove the predictions. Nonetheless, this model shows 


promise and should be pursued further. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 
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CONCLUSIONS 


1. 


The use of fins may lead to significant enhancements 
of the steam-side heat-transfer coefficient. Under 
vacuum conditions, enhancements as high as 4.0 (tube 
nunber 36) were realized, while enhancements of up to 
6.9 (tube number 37) were realized at atmospheric 
pressure. These enhancements, which were up to three 
times greater than the area ratio (finned tube 
area / smooth tube area) are due to the thinning 


effect of surface-tersion forces. 


For rectangularly-shaped fins with fin heights of 0.5 
and 1.5 mm, the 2.0 mn fin spacing appears to be the 
optimum among the fin spacings tested (1.0, 1.5, 2.0, 
and 4.0 mm). This was the case for both vacuum and 
atmospheric conditions. For Similar tubes with fin 
heights of 1.0 and 2.0 mm, Georgiadis 7 7) found the 


optimum spacing to be 1.5 mm. 


The use of fins on a stainless steel (low-thernal- 
conductivity-metal) tube degraded the heat-transfer 
performance in comparison to its corresponding smooth 
tube. In this case, the poor fin effeciency combined 
with the deleterious effects of retained condensate 
dominated over the beneficial thinning of the conden- 
sate film between fins, resulting in a net reduction 


of the heat-transfer performance. 


Fin geometry affects the extent to which the surface- 
tension forces thin the condensate film on the fin 


surfaces. Parabolically-shaped fins were shown to 
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naximize the thinning effect of surface-tension 


forces. 


Increasing steam velocity (i.e., vapor shear) from 2 
to 8 m/s increased the steam-side coefficient by 
about 10 percent (under vacuum conditions), in 
contrast to Yau et al. [1], who showed a HO percent 
increase when steam velocity was increased from 0.5 


to 1.1 m/s (at atmospheric pressure). 


For tha "optimum" tube with rectangular fins and fin 
height of 1.0 mn (tube number 6), the overall heat- 
transfer coefficient without the insert showed an 


enhancement of 40 percent over the corresponding 


smooth tube. With the insert in place, the finned 
tube showed overall enhancements of about 155 
percent. 


Most theoretical models developed thus far to predict 
the heat-transfer performance of steam condersation 
on horizontal finned tubes neglect one or more 
phenomena of major importance. Thus, there are no 
nodels currently available to accurately predict the 


performance of finned tubes used to condense stean. 


B. RECOMMENDATIONS 


5: 


Make further refinements to the Webb et al. [ 253 
model in an effort to predict the experimental data 


within a reasonable accuracy. 


Manufacture and test additional finned tubes with 
parabolic and spiral fin shapes with different dimen- 
Sions (i.e., fin spacing, height, and thickness). In 
this manner, determine the optimum dimensions for 


each fin shape. 
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Manufacture and test a series of tubes with varying 
tube diameter to study the effect of diameter on the 
heat-transfer performance. For this purpose, select 
the optimum, rectangular-shaped fins (s = 1.5 mn, 


e = 1.0 mm, and t = 1.0 mn). 


Collect data on finned tubes fitted with drainage 
strips made out of both solid and porous metal to 


determine the enhancements they provide. 


Continue tests on  finned tubes with tube-netal 
thermal conductivities between those of copper and 
stainless steel, such as copper-nickel alloys and 


aluminum. 


Conduct data runs at atmospheric pressure with steam 
velocities of 0.5, 0.7, and. 1.1 n/S, using tube 


number 6 to compare with the data of Yau et al. [1]. 


Perform data runs using different fluids with widely 
varying fluid properties (1.2% surface tension, 
viscosity, etc.) to further study these effects, and 
to compare with existing theoretical modeis and 


experimental correlations. 
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ELUCEDUEEUTUK”USE OF THE VEUBBTET XL. (257 HODEL 


As proposed by Webb et al., solving for the steam-side 


heat-transfer coefficient proceedes as follows: 


1. 


First, a value for £ must be determined. This 
parameter, as shown ty figure 2.4, is a measure of 
the fin aspect ratio (e/t), and it must be deter- 
mined using an iterative method. It is a function 
of fin geometry only, and should match the fin side 
profile as closely as possible. Webb et al. (25) 
provide details of this method, aba no further 
discussion is provided on this since a true £- value 
does not exist, as defined by Adamek [27], for a 


rectangular-shaped fin. 


With the value of known, determine the  heat- 
transfer coefficient for the fin sides in the 
unflooded portion of the tube. This value is 


obtained using the following equation `25, 27]: 


1/4 
poss 2.149 s| hfe Om Sm of (5 + 1) | (A.1) 


f “m HS r (Ë + 233 


Determine the heat-transfer coefficient for the 
unflooded root area between fins using the Nusselt 
equation written in terms of the film Reynolds 
number. For this purpose, use the following itera- 


tive nethod: 
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Calculate the condensation rate (m) 





using: 
m =n h, Ae, AT (A.2) 
hfg 
b) Assume: 
کے‎ (A .3) 
C) 


Calculate the film Reynolds number by: 


A ጣ i 
Re, x E 
f ur [ls t 5) በብ 5) 


d) Calculate h, using: 


2 -1/3 
5 Hf 
h > 1: Re 
: E 57315 (5-5) 


e) Calculate the base heat flov by: 


r he (A.6) 
f) Solve for m using: 
0 
m = + FT (2-1 
r T» ) 


g) Repeat steps c thru f until convergence is 
obtained. 
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Determination of the heat-transfer coefficient in 
the flooded part of the tube requires that a 
computer solution be used to solve the two- 
dimensional conduction problem of heat flow through 
the fins and retained condensate. This solution 
provides a ratio ($9) of actual heat flow into the 
coolant divided by the heat flow that would result 
from fins of zero thickness. The heat-transfer 


coefficient is then determined by: 


k 
LES (3.8) 
The fin efficiency is determined next using the 
following equation: 


- ۱۱۰) 1 
m (e + t/2) (A.3) 


vhere 


2 
` E (A. 10) 
Ko 


It should be kept in mind that this fin efficiency 
Poon Valid Cor fb” fins. Nonetheless, its use 
is justified since there are no equations availatle 


specifically for thick" figs. 
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The heat-transfer coefficient for the entire tube, 


then, is obtained by substituting the above quanti- 


ties into the following equation: 





= = Akt A ип ፦ 
ከ ከ 0 Ë 7 + nh E | V TTE Y 
T 


sf f As f T hz (A. 11) 
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MODIFIED WILSON METHOD 


The "modified Wilson method" calculates the leading 
coefficient for the Sieder-Tate equation. This method 
assumes a form of correlation for both the water-side and 
steam-side heat-transfer coefficients, each of which 
contains a coefficient to be determined by iteration. In 
the past, a Nusselt-type equation was used for the stean- 


side heat-transfer coefficient as shown below 77): 


3 1/3 
ED - 8 kf pf (pf y Py) һер 
с uf Do q 


(B.1) 
Equation (B.1) results in a 6 value of 0.655 for conditions 
where no vapor shear is present. The presence of vapor 


shear generally results in a higher value, which must be 


determined iteratively. The use of this equation should, 
therefore, be restricted to relatively low vapor shear 
conditions (i.e., vapor velocity less than 1.0 m/s). Note 


that a Nusselt-type equation does not represent the "correct 
form," as it contains no terms involving vapor shear. To 
avoid this deficiency, a correlation developed by Fujii and 
Honda (29) which takes the effect of vapor shear into 


account was used, and is shown below: 


175 


Nu 
= 0.96 F 
Retp : (552) 
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Pewriting ecuation (3.2) to express h as a function of heat 


flux and vapor velocity results in: 


: 1/4 1/8 

h = 5 fg 3 3 = 

o ሠ ES ۴ p Pg a .. = (8. 3) 
The constant g in equation (B.3) is to be determined by 
iteration. A Sieaer-Tate-type equation is used for the 


water-side heat-transfer coefficient, as shown below: 





i / 0.15 

i ^T s c ne Á PER i Pie. = ( 0 (8.5) 

kr Hw ` 
Equations (B. 3) and (B.4) are then substituted into the 


equation for the overall heat-transfer resistance given 


below: 











| 
+ 
十 
E 7 


(B.5) 


A linear equation used to generate the Wilson plot is 


obtained, and shown below: 


ә 
едәр " 
1 
yoo 
X ] 
Y = 2 
7ج‎ (B.8) 


The parameters X and Y are determined from fluid-property 
values and heat-flux measurements taken during data runs 
using an uninstrumented smooth tube. Iterations between the 
coefficient C, and the coefficient 8 , are continued until 
convergence of the coefficients (between two successive 
iterations) within 0.1 percent occurs. The siope of the 
Wilson plot generated is the reciprocal of the desired 
Sieder-Tate coefficient. 


APPENDIX C 
LISTING OF RAW DATA 


The following pages contain raw data obtained for tubes 
number 26 thru 44 under vacuum conditions and at atmospheric 
pressure. Raw data for tubes number 1 thru 25 are presented 


by Georgiadis [7 }. 
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UNCERTAINTY ANALYSIS 


There is always an uncertainty associated with any meas- 
urement which is dependent on the measuring-device accuracy, 
as well as on the operator's experience. Numerical data 
collected during this thesis effort were used together with 
theoretical formulations, so final values of the steam-side 
heat-transfer coefficient may be distorted due to  uncer- 
tainty propagation during calculations. In cases where the 
final results show large uncertainties, it may be unwise to 
draw any conclusions. Instead, the apparatus and/or the 
measuring techniques may have to be modified. 

Estimates of uncertainties for this investigation were 
obtained using program "UNA6" which is listed at the end of 
this appendix, along with selected uncertainty evaluations. 
This program determines the uncertainties using an equation 


proposed by Kline and McClintok [32] shown below: 


/ 
2 2 | 
Np = W + oF. U ә == ОК 1 (D.1) 
| 9:1 | 9:2 ° s 


where 











W, is the uncertainty of the desired dependent 


variable 


Xi, Xə, ==-=-, Xp are the measured (independent) 


variables 


Wi, Woe «ea, ¥, are the uncertainties in the measured 


variables 


A complete discussion covering the developement of the 
uncertainty analysis used for this investigation is given by 
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2020  Jo*Up/Lmtd ' ا۷۸۹۵(‎ SEAT TRANSFER COEF. 

2030 ጳ1።"015”።(1521”1ር0)/((1-[ር1)ታ።(15።1ር0)”።ህር((15-1ር1)/([5-1ርዕ)›)) 

777 ፐ..፣፤ !/፡.5>ሙ1ር1)5,((15-ር1)/415=1ሮዕ›)) 

2054 2R3əbtro/(( Ts-İcoleLÜÜĞ(( 1s-İci)/(Tse-Tco3)) 

2151 | dc nas R 202 72943 ) 

2079 JDuosUo*SOüR(C(XGDap/0p» 2*(Dimta/Lmta? ^2») 

2080 M=Mu 

1. 1 ኤፌ ር ኑሩ 7239 آ۹‎ 

«100  KusFNKu«T1) 

2[10 Bes). * [NTERSCEPT FROM SIEDER PROGRAM 


ہے 
Po‏ 
- 
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2.060325. ? ኒኒ ቪ ክበ ዓን Ge UNE INNES Leer IEE 


] 
2 (22.034925 ' LENGTH OF (NFINHNCO RIGPT PaWRT JF Tee 
2140 101e“.01305 
2:50” 0287750 
2160  Pr=Cou*Mw/hkw 


2170  MuysrNMuut T) 

21:80! UNCERTAINTY OF INSIDE HEAT-TRANSFER COEFF. 

003265 ST r 

2200 FHrazeoku/Dr- (CreRe 84295 4964 ነ“ ህሮ) 

2210 ወ1120/(ኮ1*።01”([+[. | ወ * 2 ዝየፍ2)።1) 

2220 Cfc-«XMuy/FNMuutc T*0t:2)^.14 

2020 IF ABS CEE 0 OL THEN 

2240 (fa((f4Cfc)e,5 

50 51 200 

260 END IF 

2270 Pi=P1-(D1+D1) 

2280 “BleXD1-01:“P 10) m5 

2290 ዛ1።(ዛ1“ዎ1/(እአእር*81)) ,5 

2200 "Pete ۷۵0 

2310 B2=(D2-Di)=PI1=(Di+D2)=* 5 

2 M2=(H1-P20/(Ke»B2))',S 

2330 بزع تمع‎ a EID AS 

2340 Fe2-FNTanh(M2-*LZ22/(M2*L2) | 
2350 О:с=9/(Р1ӱ"О1*С ет *Ее192сС2*2е2)әНн|) 

200 “7:25. (DEC De иси ۱۱٦۹9۹9 ٘ ٤+) ٰ ب01٦‎ 

2370 Dkw=.0010 * ERROR IN WATER THERMAL CONDUCTIVITY 
2280 De12.0005 ! ERROR IN SiIEDER-'ATE COEFFICIENT 

2390 Dor=.05 * ERROR IN PRRNDT. NUMBREF 

2800 Dc+=z0.E-5 

2410 ልጳ4= , 1ፅ*ዐር፣ /ር፣ 

2420 ሀካ1*ኮ1”508((ሀሄህ/ኗፍህ) 2*ቀ(ሀ0/1/01) 2+¢.3*Dre/re) ጋ2 :ዎህ፣/የያቄ፣ ሠ ደከርጅ፡:፡በ.- 
ac 

2430! UNCERTAINTY ОР OUTSIDE HEAT-"RANSrFER COEFF. 

2440  RwsDos-LOG«XDo/D12/€(29c0) ' WALL RESISTANCE 

2850 Ho CI 7loY-6Dosc/CDL «CL «Els ad Ee DD ህህ: -Ru) 
2450 DrusRu*SORCIGa20/00?^ 2* CD c/kc»"72*C0do/(Cbe-LOGCXGoo/D120)0»- mae 527 
192399825 

2470 285-1/1)0-Ru- (0o/ (0:“Hi7) 

2420 AB=Du0/ (U0 2*A5) 

2490 QB/7süru/A5 

2550 a8s«cCDo/Di1*H122*(Dhi/H1)2»/85 

2510 PRINT 

258 Dha=Ho*=SAR(AG 2+A772+98 2) 

2530' CALCULATE THE % UNCERTAINTY IN Ho 

2540 Ргһо=9Һо" (00/Но 

2550? CALCULATE TRE % UNCERTAINTY IN REYNOLDS NUMBER 
2550 Prre=Dre*100/Re 

2270! CALCULATE TRE 4 UNCERTAINTY IN MASS FLOW RATE 
2530 Prmt=Dm**=100/Mf 

2530' CALCULATE THe 4 UNCERTAINTY IN HEAT TRANSFER 

2600 Prap=Dap*100/Up I 

2610! CALCULATE THE % UNCERTAINTY IN LMTD 

2620 Primtd=Dimtd-100/Lmtc 

2630! CALCULATE THE 4 UNCERTAINTY IN Rw 

2640 Prru=Dru*100/Ru 

2650! CALCULATE THE Z UNCERTAINTY IN OVERALL PEAT TRANSFER CODEF. 
660  FruesDuos100/Uo 

26/70! CALCULATE THE % UNCERTAINTY IN INSIDE HEAT TRANSFER COEFF. 
2680 PrhrsÜhi*1007/H: 

290 PRINT 

2700 PRINT USEING “T0OX.°°UNCER TAIN, Yee Nate SS ርር 
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PRINT 

POINT USING IU," VARIABLE 

PRINT 

INT USING IOUS. “ass Flow Rate. Ma 
PRINT USING "10X,""Reynoias Numoer, Re 
EwIMI USIMG "10X.""Heat Flux. 


q 
157 07” oq”"ean-lon ÜL f. LMİD 


PRINT USING “10X. "Wall Resıscance, Ru 
PRINT USING "10X.""Üverall H.İ.Ü.. Uo 


SLN USING “TUX. “Water-srde Ҹ.1.С.. Hi 
PRIN US NG VOX. Steam-Side H.1.C.. do 


END 

DEF FNMuu(CT?) 
82247.8/€7*133.15) 
Muus2.4E-5«410^A 
RETURN Muu 
FNEND 

DEF "NTann(X)? 
Р-ЕХР(Х) 
(-2ХР(-Х) 
тапһ=(Р-О) /(Р-ј)) 
RETURN Tann 
FHEND 

DEF FNKu( Tt) 


FESCENTTONCSS DENIS T 


7500: am” 
КР О” ВБутв 
715 r 
ә ی٤‎ 
BOIL. PER 
۰00 200 ruo 
Да ри ири 
DE oO. Te TNO 


25 . coco) 194). S00/-11=(.52577- .073944"=)1))) 


RETURN hu 

-MEND 

DEF *NMutCT) 
Q«247,3/€T241353. 15)» 
4y=2,4E-5-10 4 
RETURN Mu 

ЕМЕНО 

DEF 2ҸЕҺосТ) 


5:526:177:55::- ES 0 ə CE 3 ii. 2416785?) 


RETURN Ano 
MENO 
DEF FrNCpu(CT) 


бидә 1106 1442 6026Е-2-1*(43.423612-54*2,71428Е-7әТ)))*1000 


RETURN Cpu 
EHEND 

DE* rHTvsv(CEmt» 
እ በከን OT? 
T=C (0) 

FOR 1=! TÜ ? 
Т=13С(1)"ЕлР 1 
НЕХТ." 

RETURN T 

FNEND 
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DATA FÜR THE UNCERTAINTY ANALYSIS: 

“6 IE racueb 
Pressure Cardi ۹۰۶ AON TLE) 
ጋቱ MEETS - 2 ٦) 
Water 550657: = 79,00 
Water Velocity = NI Oum, 
Heated = QUNN COS (H/m 2D 


Mass Flow Rate ٣ 
Revnolas Number, Re 
Heat 0۳ C 
1510 M I 
all EEE ETE ин 
ara RM Q 

:7 تمت رك -ع تخ جيم 

“raam A 


dE R E 


Pine ESTE 
امم مت عتم م‎ Лора 
Steam emas ra mu = 
Wa ә, Flow 77 
mater Yeioclty = 
feet O = 
LUNCERTAINTY ANALYETE: 
VARIABLE 
Hass 0006 
Reyno. de Number. Re 
Heat Fiux. q 
Lgs amnes mM O 
ual? Kola OL md 
Hera ጸን. )) 


Water-Side HM. 
“tean. tae 
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ANALYYSI 


CE cnm 


ERAL 
2. 
тб 
122 
2.57 
127: 
4.46 
—. c» 
و کے سے‎ m 
راد‎ 
52 اناا‎ 007 ES 
28.79 (eg () 
3 
3.31 (nas) 
3.808£«05 (l/m 2) 


E ٠٣‏ ۲م 


CD 一 一‏ رح 
C‏ 
= 


[3 መሙ 
OO Ls 


DATA Ја HE UNCERTAINTY SNALYSIS: 

Fiie tham e: r 3255 

ከሙ ጉጥ “ሥሎ. wond i Pen: Atmospheric “101 ርን 
2755" iemperature xu, eau) 

NST Flou Rate (%4) = Q20.00n 

| ከከ ከ Velocity iim la/s) 

255 ”lux = £,502E*05 (H/n 2) 
UNCERTAINTY ANALYSIS: 

VAR TABLE OCF RCEN™ UNCERTAINTY 
fess r Qu Rate. Md =. 45 
Reynoids ПНипбег, R 550 
Heat = lux. a m 
кә ean -Tam Duc. 0 ^. 07 
Hali Resi: stance, ጂኒ) 2.57 
nari n. TEN. llo 1 eu 
keter Side 205. Hi 2 
ےم مر دیدج ہت‎ FDİ.U,., Ho 26.08 
DATA OR THE UNCERTAINTY Ant YS r$: 

:7.::—-: 2َ ከ 
E erc ء-ے-‎ 3 ron: Hun.” a) 
ከ.  “5ertarture 775: (deg C) 

uer Eis tere (7) = б0.00 

Water Veiccity = zr Ac) 
277765 = |.164606 (H/m 23 

£1 5:07 ጋ በ i 7215: 

VARTABLE OCRUDMT UNHCLRTATNIT 
Маз= Flou Kate, Md 1 ,44 
Reynolas tlumber, Ke fees 
016 iux. e 0.48 
570 61 Dit. LATD 0.14 
us ede, Ku 2.57 
eran | H.T, 07 To Си 
Hater de HO... Hi 1101] 

Steam- Sige H.T. 300 HO 2.80 
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DATA FOR (HE "+۵7 


т‏ ا 

res sure Condition: 
Stean enpera S 
Water Flou Rate (7) 
Water Velocity 

Heat 210 


۸ 1 ۶٣ 
VARIABLE 


Маса lom Eu a 
reynaids Number, Re 
A EE 
ou GI Ül 
Wal! Kesistance. 3 
Пуер ур 2 ፡ (የሙ. 
Water-Side H 6 

S vede 169 C DE W 


392 


2 


US cuun (11 Pap 


48 


50 (Дев СӘ 


¿O Om 


4 


1 


5 


Lo (my => 
73E*05 (W/m 2) 


PERLENT ۱۹ ۶ 7 


DATA POR TBE ٣ sS 


ки ии 

ወሮ. Cae ከ0: 

“ ٣ 3۳۰ یھ‎ ٣ = 
water Flou kale m) = 
Water Velocity = 
Heat Fluy = 


n 
(c 
ад 
2Ј 
ad 


ተ 7 
VARIABLE 


Masa Flou ةا‎ 
Raynoids Number, Re 
8 3٠. 
|оа-Неап-Теп Diff, LHTD 
Wali Resistance, Ru 
Overa Mn o 
Water Sio- CS ይ. 
Steam cide ии 
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1 


` ә уран ин 
ራሩ Cd کے‎ 02 O 


“一 =} 


= 


D 


Uf) 


I 


—2 1 Flu fi 


- 
` æ 


Pope پ پھل‎ 
(905 €$C3FP9 1 
DO CYL DEIS 


C 


7” 
Uum Cd 1050 
¿41 (eg ©) 

BO 

55002 ٦ 
۵)٢) ከ 1... 


PERCENT UNCERTAINTY 


7- 
516 
م‎ n) 


-30 
45 
39 


Is о O cci 


0١۱ ز٠٦‎ . 
D N Ù 


popu n HE USCERTInIMIY ANALYSIS: 


File flame: £ ۶۴ 

“Pressure Condition: Atmospheric (10° kPa) 
oteam lemperature = Cec C) 

Newer Flow Rate (7) = 20.00 

Water Velocity = oec) 

Heat rlux 2577 : 05 (MR 
UNCERTAINTY ANALYSIS: 

JAR TABLE PCRCENT UNCERTAINTY 
Mase Rate. Md 
Reynoids Nummer, Re 5: 

Heal F Ure -g um 
ES ean lem Dit. MID Huk 
ASS is tante Ki 7: 
ከሽ ን ከ. 6+ Uo aca! 
ር 51-5108 myleC., me m 
ይ ə 


ከ ስ”! ل‎ Ho :6.b? 


ВЕ ۱٢۷۷ ۹۹۳31۲۲۴ ۱۸۷5 ۱ 

pore ህክ. 

Pres Mre EIS Ton: Rinotpneric (im s) 
Stean .emperature = 49,79 (Dea ©) 

T E (23 = Gt 07 

Water Velocity = 77 4010010) 

ር = 1,409 +ሀፍ (ህጩ/ጦጠ 2) 


AE ERT AINI GHAEYSIS: 


VARIABLE CERCEHT. 'INCECRTATNTY 
15-00 1000 pare, Md 1.42 
Reynolds tumoer, Re 1.24 
nest 7707774 0.48 
Loa-“Mean-l=m Diti, LIMID kə s 
Wall Resistance, Ru 457 
کے‎ a l mm mtns 59 . 
Water በር |. Hi 1.46 
STEAM IIE H.T.L.. Ho 575: 
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